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SUMMARY 



Thermal-electric means of preventing ice formations on 
a propeller "blade have "been investigated and a theoretical 
"basis for the continued development of thermal-electric 
blade shoes is provided, 

A method is presented that can "be applied to the design 
of thermal-electric "blade shoes for any propeller or rotor, 
and c?.n optimum heat distrihution is determined for a propel- 
ler "blade. 



INTRODUCTION 



The design and development of suitable thermal ice- 
prevention equipment for aircraft propellers have been under- 
taken as a part of a general research program concerning ico- 
provention equipment for aircraft o 

The National Research Council cf Canada has conducted 
flight tests, under natural icing conditions, of propellers 
equipped with electrically heated blade shoes (reference l). 
These results formed the basis of the preliminary blade-shoe 
designs used during flight tests conducted by the NACA in 
the vicinity of Minneapolis, Minn*, during the winter of 1942- 
43 (reference 2). The blade shoes tested by the National 
Research Council of Canada and the NACA consisted of a layer 
of neoprone serving as thermal and electric insulator bonded 
to the propeller blade with an outer layer of electrically 
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conductive neoprane tonded to the insulating layer. (See 
fig. 1.) The power was sux)-Dlied to these blade shoes throu,-^h 
radial vrires "between the neoprene layers at the rear edges 
of the shoe. Flight tests with electrically heated propeller- 
blade shoes reported in references 1 and 2 have established 
the practicability of protecting propeller blades from ico 
formations by thormal-electric means. The data of references 
1 and 2 indicate that the power required for ice prevention 
may be excessive for certain applications, although suffi- 
cient pov;er for some degree of ice removal may be provided 
readily. The analysis reported herein provides a rational 
basis for establishing thermal-electric pr opeller-blade- shoe 
de si gns . 



SYMBOLS 



The follov/ing symbols have been used in this analysis: 
A aspect ratio 
B number of propeller blades 
D propeller diameter, feet 

J work equivalent, 778 foot-pounds per Btu 

ITug boundary-layer ITusselt number based on the laminar 
boundary-layer thickness, h6 /k 

P normal pressure coefficient on airfoil section 

normal pressure coefficient due to the additional lift 
di st r ibut i on 

P-b normal pressure coefficient due to the basic lift 
di str ibut i on 

Pp reference profile pressure coefficient 

Q, heat quantity, Btu per hour or watts 

R propeller radius, feet 

Reynolds number based on blade chord, ^^^p^c/v 
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Reynolds nurnlDer "based on leading-edge diameter, '^h^/'^ 

rate-of- icing ratio 
S area, square feet (unless otherwise noted) 
T temperature, degrees Fahrenheit 

Tj^ t Giiiperature rise due to aerodynamic heating, degrees 
Fahrenheit 

V velocity just outside of the boundary layer, feet per 
sec ond 

Vj^ maximum local velocity, feet per second 

local velocity due to airfoil shape, feet per second 

Vq forv/ard velocity of propeller, feet per second 

resultant velocity of the blade section, feet per second 
local velocity on the reference profile, feet per second 



local velocity at the laminar separation point, feet 
per second 

W- water content of the air stream, pounds of water per 
cubic footofair 

a slope of the lift curve for any blade section 

slope of the lift curve at infinite aspect ratio 

"b chordwise distance through which heat is conducted, feet 

c blade chord, feet 

Cj propeller-blade section lift coefficient 
cj additional section lift coefficient 



ib 



basic section lift coefficient 



Cp specific heat of air, Btu per pound, degrees Fahrenheit 
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Cp specific heat of water^ Btu per pound, degrees Fahrenheit 
d leading-edge diameter, feet 

g acceleration due to gravity, 32.2 feet per second per 
sec ond 

h heat-transfer coefficient, Btu per hour, square foot, 
degrees Fahrenheit 

k thermal conductivity, Btu per hour, square foot, degrees 
Fahrenheit per foot 

I distance from the edge of an ice layer, feet 

m distance normal to blade-shoe surface, feet 

n propeller speed, revolutions per second 

p absolute value of the ratio of the slopes of a double- 
roof velocity profile 

q heat per unit area, Btu per hour, square foot or watts 
per square inch (as noted) 

r blade station radius, feet 

8 chordv.'ise distance along the airfoil surface from the 
stagnation point, feet 

t thickness, feet 

X distance along the airfoil chord line, feet 
$ helix angle, degrees 

helix angle at 0.75R at (^o/^^^max' cLegrees 
Q angular velocity of the blade, radians per second 
a angle of attack, degrees 

ai angle of zero lift, degrees 
*' 0 

3 blade angle measured from the plane of rotation, degrees 
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8 laminar "boundary-layer thickness from the surfac 
the point at which V = 0.70? Vp^, feet 

8rp heat-transfer characteristic length for a turbul 
boundary layer 



i 


turbulent boundary- layer pararaet er 




e 


blade anvq:le at any radius, de(^;rees 




6 m 


blade angle at 0.75 R at ("^o /^^ )max » 


degr ee s 




kinematic viscosity, square feet per 


se c ond 




ice-thickness parameter 




Subscript s 




I 


insulat i on 




L 


lower surface of the blade 




U 


upper surface of the blade 




\) 


blade 




h 


blade root 




i 


i ce 




s 


laminar- separat ion point 




t 


blade tip 




u 


leading edge of the blade shoe 




V 


rear edge of the blade-shoe conducting layer 


w 


trailing edge of the blade 




a s 


lolade surface aft of the blade shoe 






blade shoe 




cl 


conducting layer 




c s 


blade-element radial cross section 
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stag stagnation-pressure region 



0 



am"bient air 



1 



a point on the "blade or shoe surface 



2 



a point in the "blade shoe 



0-75 r/H = 



Oo 75 



+ 



positive values of I 



negative values of I 



METHOD 



An optimum design of a propeller-blade shoe should pro- 
vide equal protection to all points on the blade-shoe surface 
at the design conditions and have the minimum permissible 
total power input. The method of analysis presented herein- 
after is general in scope and can readily be applied in the 
design of an optimum blade-shoe arrangement for any propeller 
or rotor blade. 

In this analysis the blade shoe is considered to extend 
radially from the -or opeller-blade root, station (r/H)j^, to 

the T)r opcller-blado tip, station (r/E)^, and chordwise from 

the blade leading edge, station (s/c)^, to station (s/c)^ 

on both the upper and lov/er surfaces of the forv/ard portion 
of the blade. (See fig. 1.) Heat is considered to be ap- 
plied to this portion of the blade while the after portion, 
extending from station (s/c),^ to the propeller trailing edge, 
station (s/c)t.;, is neither heated nor covered. The blade 
shoe is considered to be comprised of an inner layer, adjacent 
to the propeller-blade surface, of an insulating rubberlilie 
material of thickness tj overlaid with an electrically con- 
ductive layer of similar material of thickness t^. The 
use of an additional thin layer of insulating material over 
the conducting layer to increase the resistance to abrasion 
will not be considered in the analysis; however, the effect 
of such a layer can be determined readily. 

Ice protection may be accomplished by two different 
processes: (l) ice can be prevented from forming, and (2) 
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ice can "be removed periodically in thin layers. To provide 
for ice prevention on propeller "blades, it is necessary to 
supply sufficient heat to maintain the temperature of the 
hlade-shoe surface ahove 32^ F in any atmospheric condition* 
Protection "by ice removal can "be effected iDy providing suffi- 
cient heat to raise the "blade-shoe surface temperature auove 
32^ P only after ice has formed. It can "be seen that the 
removal process involves thermal economies by utilizing: the 
insulating properties of ice to reduce the convection losses 
and to ohviate the need of supplying the heat of fusion and 
the heat of evaporation. 

The total power required for either ice prevention or 
ice removal is dependent upon the atmospheric conditionsi the 
propeller-operating conditions, and the propeller design, and 
may "be expressed as the summation of the heat required to 
prevent or remove ice formations on the propeller-blade-shoe 
surface and the heat lost through the after or uncovered por- 
tion of the propeller "blade. The heat req^iired to prevent or 
remove ice on the "blade-shoe surface Qv^ g is the summation, 

over the radial and chordwise extent of the "blade-shoe sur- 
face, of the values of unit hoat (heat per unit area) re- 
quired for the prevention or removal of ice at each point on 
the blade-shoe surface q'bs* Lil^ewise the heat lost through 
the after surface of the propeller blade 0,^3 is the summa- 
tion, over the radial and chordwise extent of the tmcovered 
portion of the propeller blade, of the values of unit heat 
lost at each point on the after portion of the blade Q.as* 

The unit heat required to prevent or remove ice at each 
point on the blade-shoe surface qi^g is dependent iipon (l) 

the surface neat-transfer coefficient, (2) the water content 
of the air stream, (3) the ambient-air temperature, and (4) 
the resultant velocity and effective angle of attack of the 

propeller-blade section. The unit heat loss o^^ at each 

■as 

point on the after surface of the propeller blade is a func- 
tion of (l) the surface licat-transf or coefficient, (2) the 
ambient-air t emperattir e , (s) the thermal conductivities of 
the propeller-blade and blade-shoe materials, and (4) the re- 
sultant velocity and effective angle of attack of the 
propeller-blade s ect i on. 

The method of analysis which follows is therefore con- 
cerned v/ith the determination of these various factors v/hich 
are necessary to evaluate q. ^ and q and thereby effect 
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a solution for the total power required and for the optimuxn 
hoat distribution for the specific application. 



Determination of the Surface Heat-Transfer Coefficients 

The sulDject of heat transfer from airfoils has "been 
widely investigated and several approaches to the suhjcct 
have "been developed. Since each section of a propeller "blade 
is essentially an airfoil section, the method of determining 
the surface heat-transfer coefficients of an airfoil section 
has heen utilized in the solution of the surface heat-trans- 
fer coefficients of a propeller "blade. 

The surface heat -trans fer coefficients of an airfoil 
section are dependent upon the physical characteristics of 
the airfoil "boundary layer, and therefore its determination 
requires a knov/led^^ie of the location along the airfoil sur- 
face of the point of transition of the "boundary layer from 
laminar to turoulent flow as well as of the "boundary-layer 
thickness along the airfoil surfs.ce, 

P r opeller-hlade section velocity di strihution , - Since 
the extent of the laminar boundary layer for a specific air- 
foil section - or, as in this analysis, a specific propeller- 
hlade section - is dependent upon section velocity distribu- 
tion, one of the first steps in the determination of the 
propeller-blade surface heat-transfer coefficients is the 
determination of the velocity distribution for the propeller- 
blade sections under consideration* The velocity distribu- 
tion for any propeller-blade section is determined by the 
lift coefficient and angle of attack at which that blade 
section is operating^ The angle of attack of each propeller- 
blade section can be v/ritten as the difference betv/een the 
blade angle and the helix angle at the blade station consid- 
ered , 

a = 0 - $ (l) 

and the lift coefficients for each propeller-blade section 
can bo determined by correcting the section lift-coefficient 
curves for the effect of finite aspect ratio. The effective 
aspect ratio for each blade section can be obtained from the 
foil owing equat ion: 

. (2) 
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which is developed in reference 3, and in which 

Be 



^0 = 



ttE 



From the effective aspect ratio, the correct lift-curve 
slopes can be computed liy the equation 



a = 



1 + 



a^ X 57.3 



ttA 



(3) 



and the correct lift coefficients obtained hy the relation 

cx=a(a-ai^) (4) 

If the lift coefficient and angle of attack are known, the 
velocity distribution over each pr opeller-hlade section can 
he determined by the method of either reference 4 or 5. 

Laminar bound ary layers . - The extent of the laminar 
boundary layer is a function of the blade-section velocity 
distribution and the boundary-layer Reynolds number. For 
the purpose of this report, the boundary layer will be con- 
sidered laminar from stagnation to the laminar- separat i on 
point and then turbulent to the trailing edge. The method 
of locating the laminar- separat i on point is developed from 
the theory by Von Karman and Millikan (reference 6). By 
expressing the relation of reference 6 in the nomenclature 
of this report, the following equation is obtained: 

2 

(5) 

in which X is equal to Gg of reference 6 and is explic- 
itly defined therein. The absolute value of the ratio of 
the slopes of the double-roof velocity profiles, which are 
drawn to approximate the actual velocity profiles, is p 
and the curve for Vg/Vjvj as a function of p is plotted in 
figure 2. 



M 



1 - p(X - 1) 
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The heat transfer in a laminar "boundary layer is shown 
in reference 7 as "being dependent on the "boundary-layer 
Nusselt nuin"ber and is a function of the "boundary-layer thick- 
ness 5 and the thermal conductivity of air k. In the re- 
lation 



Nug is a function of the shape of the "boundary-layer 
velocity profile and Prandtl numher and has a value of 0.765 
for the 31asius-type "boundary-layer velocity profile for air. 
Experiment has shown that, over the forward portion of an 
airfoil where a favorable pressure gradient exists, the 
Blar.ius-type h oundar y-layer velocity profile is closely ap- 
proximated. At points on the surface of a "body downstream 
of the minimum-pressure point, a laminar "boundary layer ex- 
hi"bits a tendency to separate. Since the velocity gradient 
at the surface decreases as separation develops, the value 
of ITug must diminish until at the separation point its 
value is zero. It is cnn^ldered that reducing Nug linearly 
from the mi nimum-pr e ^5 r e point to the separation point will 
satisfactorily appr o 7.imat e the actual cas'2. Accordingly, in 
the anal,ysis of this report Hug has leen considered to have 
a constant value of 0.765 from the stagnation-pressure region 
to the minimum-pressure point (maximum velocity) and then to 
vary linearly to zero at the laminar separation point. 

The laminar "boundary-layer thickness at any point s^, 
as used in reference 7, is given hy the equation 




(6) 




(7) 



8o 17 




C 



which can he changed to the form 



6 = 




0 
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to simplify the computations. Equation (8) is theoretically 
applicaole only to Blasius-type boundary-layer profiles 
(stagnation to minimum-pressure point), "but, as noted in ref- 
erence 7, experim.ent has shown that little error is caused 
"by its use to the lami nar- separ at i on point. 

With the ahove method the solution for the boundary- 
layer thickness at the stagnation region "becomes indetermi- 
nate; therefore, another method must "be used to determine 
the heat--t ransf er coefficient in this rei?:ion. The expression 

6 ^ = c / J-.^i^ (9) 
stag J ^c.y 

c 



also from reference 7, gives a satisfactory value for the 
"boundary-layer thickness in this region "based on the local 
radius of curvature F, v/hich is not necessarily the lead- 
ing-edge radius. 



Tur'bulent "boundary layer 



transfer coefficients for 
"be obtained by the method 



The values of the heat- 

the turbulent boundary layer may 
of reference 8. The surface heat- 



transfer 

e Quat i on 



coef f i cient 



is defined in reference 8 by the 



h = 0.76 



6t 



(10) 



and the value of 6rp, the function of the turbulent boundary- 
layer thickness, is given by the equation 

6t = \^ , (11) 



in which t is a turbulent boundar y~layer parameter. The 
value of ^ at any point on the airfoil surface dov/nstream 
from the laminar- sepatat i on point is obtained by a step-by- 
step solution of the equation 

^ 6JL3 dV ^ V ^ (12) 
dx V dx i> 
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from tlie point of laminar separation. The value of t, at 
larainar separation is o'btained from 



r 0.289 8 1 
[4.057 Vg J 



i = 2.557 lege 
Squat ion (12) has "been rewritten in the form 



r 6.13 d(v/yR) -Vec ^ 
y u" (V/Vr) d(x/c) V 



(13) 



(14) 



for ease of computation. The values of f(t) related to 

t are plotted in figure 3. A simplified method of applying 
the turbulent "boundary-layer aquations to obtain the value 
of St at any point directly is presented in reference 9. 



Heat Required for Ice Prevention 

V/hen the surface of a blade shoe covering the leading- 
edge portion of a propeller blade is heated to 32° F by 
applying povrer to the electric conducting layer, there v.^ill 
be some heat loss through the insulating layer and blade 
material to the air stream over the uncovered aft portion of 
the blade » The total pov.^er required to protect the propeller 
blade will then be that required to maintain a 32^ F blade- 
shoe surface temperature plus the heat lost through the aft 
portion of the blade surface. 

Heat r equired to maintain the blad e -shoe surf ace _at 
3 2° F . - The heat necessary to maintain the blade surface at 
32° F in icing conditions is a function of the surface heat- 
transfer coefficient h, ambient-air temperature Tq, tem- 
perature rise due to aerodynamic heating Tj^, weight of 
water in each unit volume of air W, and the resultant veloc- 
ity of each element of the propeller blade . 

The unit heat required is given by the relation 
= h(32-To-T^) + V; Cp ( 3 2-T ^ ) ( 36 00 ) - .Ol- (3600) (15) 

in which the first term is the heat transferred to the air 
flowing over the blade, if the blade is maintained at 32° 
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and the values of h are those that were determined in the 
previous section. The temperature difference "between the 
"blade surface and the amhient-air stream will be decreased 
by the amount of the temperature rise Tj^ caused "by air 
flowin^^ over the blade or being stopped in the stagnation 
region. The second term of equation (l5) is the heat re- 
quired to raise the temToerature of all the water that strikes 
the blade to 32^ 1?. The third terui is the hoat equivalent of 
the kinetic energy of the water striking the blade. The value 
of 3600, by v/hich the second and third factors in equation 
(15) are multiplied, is to obtain dimensional consistency so 
that the units of q-^g are Btu per hour, square foot. 

Studies have been made to find the relation between the 
path and size of water drops as they approach and contact an 
object in an air stream, and it has been found that as the 
air-stream velocity and water-drop size increase, the path 
deflection decreases. Since the resultant velocities of the 
propeller-blade sections r.re relatively high and since large 
water-drop sizes are associated with severe icing conditions, 
this analysis assumes the cieflection of the water drops to 
be negligible. The v;e-*.ght rate of water striking the blade 
at any point, or the icing rate, therefore, is a function of 
the slope of tne surface relative to the water-drop path at 
that point. The rate of r.cing at any point on the blade can 
be expressed as a function of the icing rate at the stagna- 
tion region.. This furction is proportional to the slope of 
the blade-section surface r..nd is termed The values of 

Rl are measured fro:.Ti the soot ion layouts as shown in figure 
1. For sections with flat lower surfaces, for points 

on the lower surface is taken to be the sine of the angle of 
attack measured from t}:e angle of zero lift. Per the upper 
surface, the error duo to measuring the slopes from the chord 
line rather than the relative wind is small and conservative 
and tends to allow for small changes in co . Sach term in 
equation (15) that expresses the effect of water content must 
therefore contain the factor to correct for the slope of 

the blade- sect i on surface. 

The temperature rise at the stagnation point on a 
propeller-blade section is due to the pressure rise at that 
point and is 



T. . = -± (16) 

2Jgc^^ 



-^stag 
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as noted in reference 10. For any point other than the stag- 
nation point, the temperature rise v/ill "be taken to "be 0.8 
of that at stagnation pressure, or 



Ta = 0.8 T 



-^stag 



(17) 



The 0.8 factor includes the effects of the laminar and tur- 
bulent "boundary layers and compressibility, as defined in 
reft^rence 8, and was approximated from the data of reference 
11. 



Equation (15), with the effects of aerodynamic heating 
and blade- surf ace slope included, becomes at the stagnation- 
pressure region 



^b( 



= h 



stag 



stag 



2 -1 



32 - - 



2Jgc 



Pa' 



(32 - T„) - 



2Jg 



(3600) 



(18) 



and at any other point on the "blade shoe 

0«8 Vr" 



32 - T - 

0 



2J6C 



Pa 



+ R. V/ 
1 



% ^r(^2 - ^) 



'R 



(3600) 



2JgJ 
(19) 



The total heat required at the blade-shoe surface to 
prevent ice is obtained oy the following integration which 
can best be performed graphically, 



bs 



= BR 



(r/H) 



(r/R)h 



V 



=u 



d(^-)(20) 



where (^/^^)h» (^/^)t» ^u » c onduct i ng-layer 

limits of the radial and chordwise extent of the blade shoe. 



Heat los s.- In order to provide sufficient heat to main- 
tain the blade-shoe surface at the denign temperature, an ex»- 
cess of heat must be supplied to the blade shoe. This excess 
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of heat represents the heat lost to the atmosphere through 
the thermal circuit comprised of the blade-shoe insulating 
layer, the "blade material, and t he tiade- se c t i on boundary 
layer. At any point on the after surface of the propeller 
blade, the heat loss per unit area qg^g will be a function 
of the over-all thermal conductivity of the heat-loss circuit 
(k/t), the temperature at the point of maximum thermal poten- 
tial T3, the ambient-air temperature Tq, and the tempera- 
ture rise due to aerodynamic heating Tj^, and may be ex- 
pressed as 



The location of the point of maximum thermal potential 
in the conducting layer of the blade shoe and the temperature 
at that point will be determined by the relative magni- 

tude of the heat transmitted to the blade-shoo surface and 
the heat lost to the blade. Since is the maximum tem- 

perature which exists in the blade-shoe conducting layer, a 
decreasing temperature gradient will occur between the loca- 
tion of this point of maximum temperature Tg and the sur- 
face of the insulating layer. The distance through which 
this decreasing temperature gradient exists in the conduct- 
ing layer can be considered as an increase in the thickness 
of the insulating layer. This effective increase in the 
insulating-layer thickness will be termed Atj and is illus- 
trated in figure 1. By assuming the temperature gradients 
from the point of maximum thermal potential to the blade-shoo 
surface and to the surface of the insulating layer to de- 
crease linearly with distance, the value of Atj can be 
expressed as 



(32) 



The value of can then te determined with sufficient 

accuracy hy the equation 

^3 = 32 + U^?] I^^l-lili) (23) 

where [^I^^qJ^^ is the average value of the unit heat re- 
quired at the blade-shoe surface at the particular blade 
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station as determined from the surface heating requirements. 

The over-all thermal conductivity of the heat-].oss 
circuit ,(k/t) can now be defined as a function of its com- 
ponent thermal resistances, the insulating-layer resistance 
(tj + Ati)/ki, the blade-material resistance h/k-^ (where 
"b is the chordwise distance through the blade material 
which the heat is considered to travel), and the boundary- 
layer resistance l/h, or 

£ = f ^ „^ ^ (24) 

t L kj k^ hJ 

The determination of the thermal loss through a blade 
element of unit width (fig. 1(a)) can be simplified by assum- 
ing certain conservative changes in the heat-loss circuit. 
The blade-shoe area and the area of the after surface of the 
blade will be considered to be twi?,e their projected areas 
or S"bs and S^g, respectively, while S-^ is the projected 
blade area. The blade- ?leinent radial cr o s s- se c t i on area S^s 
is the area of the blade material at any point Si through 

which the heat loss is considered to be conducted and is also 
shown in figure l(a). Por the purpose of the calculations, 
the value of S^^ can be used more easily by assuming the 

blade-el emont m^'.torial to be a wedge v/ith the base height 
equal to the maximum blade sccticn thickness tjj^^:: 
blade station, the base width ogtial to unity, and the length 
equal to the chnrdwi3e distarce from the trailing edge of 
the conducting layer tc the blade trailing edge, as illus- 
trated in figure l(b). The blade material under the blade 
shoe is neglected in this heat-loss approximation. The re- 
sult of these assumptions is an equivalent thermal-loss 
circuit for each point si on the after surface. As shown 

in figure 1(c), the heat loss travels from the point of 
maximum temperature T^ through the effective increase in 

the insulating-layer thickness Ati and then through the 
insulating-layer thickness tj to the blade material. Since 
the material beneath the blade shoe has been neglected, the 
circuit is represented by considering the blade-shoe layers 
to be applied to the forward face of the block of blade ma- 
terial of length b, unit width, and thickness t^g as 

given by the wedge approximation for any point s^. After 
the heat has been conducted through the block of blade 
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material, it is lost to the ambient-air stream by forced 
convection through the boundary layer. The assumed equiva- 
lent heat-loss circuit is compLete if the air stream is con- 
sidered to flow over the opposite end of* the block from that 
which is considered covered by the blade-shoe layers. The 
complete equivalent heat-loss circuit as described is shown 
in figure 1(c). 

If the unit heat distribution w^is constant over the 
blade-shoe surface, the total heat would be given by the 
equation 



But 



t/k 



^a s 



^as 



f; "as z ~ ^ 0 " -^A 



/ti+Ati 1 \ 




(25) 



'as 



where the unite of the t<-rins are the reciprocal of Btu per 
hour, degrees Fahrenheit. The purpose of tnis step is to 
express the three factors of the tiiermal resistance in their 
proper relation as determined fcy the different heat-transf sr 
areas for each term. In order that the thermal conductivity 
k/ t can he used for tha case of nonuniform heat distribu- 
tion, or in the expression for q^^ (equation (21)), the 
foregoing equation is rearranged as follows: 



k/t = 



L V 



a s 



OS 



/ b 



as 



r hJ 



(26) 



therefore, for any point on the blade after surface the unit 
heat loss can be dotcrmined by the equation 



■■as 



+ At 



2 0 



'as 



<27) 



The total heat 
the values of 
pr opellor 



^as 
blade as 



^oss q^^^ is obtained by integrating 
over the upper and lower surface of the 
f 0 1 1 0 v; s i 
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iir/R) (38) 



It has been shown that the evaluation of equation (27) 
depends on two unknowiis Atj and which are related. 

Since Atj depends on the relative magnitudes of the heat 
loss Q,a^g and the heat supplied to the blade shoe Qbs» 
initial assumption of Atj must be made and checked by solv- 
ing for Q^g and the resulting value of Atj used in a 
second solution for Qg^g until the values of Atj, Qas* 
Qt)g are in reasonable agreement. 

Heat Required for Ice Removal 

When the quantity of heat supplied to a propeller blade 
is reduced below the quantity necessary for ice prevention 
and ice is allowed to form on the propeller blades, the ice 
will be throv/n off after reaching a certain thickness, which 
is a function of the power input and power distribution. 
The insulating quality of the ice allows the blade-shoe sur- 
face temperature to rise above 32^ F and a water interface 
to form betv/een the ice and the shoe. The reduction of the 
adhesion combined with the centrifugal force on the ice due 
to its nass will result in ice removal. 

The process of ice removal can be made to occur at dif- 
ferent points and at different intervals, depending on the 
heat distribution and the local rates of icing and heat 
transfer. The ice thickness at a point on a blade station 
can be expressed as a function of the maximum ice thickness, 
which occurs in the stagnation region, and the rat e-of -i ci ng 
ratio 

If ice is assumed to accumulate on the blade shoe in a 
smooth continuous layer with no free edges and with the 
thickness distribution as given by equation (29), the unit 
heat required for ice removal will be given by the equation 
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'lbs = 



- T. - T, 



W V 



1 

h 



R 
2Jg 



(3600) 



(30) 



The foregoing- equation is analo^rous to equation (l5) except 
that the effect of the ice layer has "been included in the 

1' 



thermal resistance 



I -f 

equation (15) , the heat required 
of the air-stream v/ater content 

(50) 



and that the second terra of 



to increase 
to above 32^ 
will have a 



the temperature 
F, has "been 
value of 32^ F 



omitted. In equation vws^, 
only for the condition in which a continuous ice layer exists 
The second term of equation (30) is the kinetic heating due 
to the water drops striking the "blade and it can he neglected 
because its omission introduces a small and conservative er- 

neglect ed , • the heat required as indi- 



r or , 
cat ed 



If 



the t erm is 



^bs 



- T - T, 



1^ 1^ 



(31) 



Lk 



h. 



will not be "a function of the water content of the air stream 



In natural icing conditions ice does not form in contin- 
uous layers but rather in patches, ridges, or particles. At 
the edge of a discontinuous ice layer , there will be a heat 
less from the iced to tho uniced portion of the blade shoe 
with a temperature gradient along the shoo surface normal to 
the edge of the ice. (See fig. 4.) If a semi-infinite ice 
layer is assumed to extend from a point on the blade-shoe 
surface at which the surface temperature is T^^, the sur- 
face temxierature T, 



at any point 



a distance 
the surface 



I from the 
t emrer atur e 



beneath 
' T 



- 0 ' 

the ice 



T. 



= 1 - 



location or i q as a 
at I equals infinity 

h 



layer at 
function of 



(1 + ^o/Ti^)(l 
1 + \1; 



) 



e 



(32) 



wher e 
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^ = 



k. 



(33) 



which is developed in the appendix and illustrated in figure 
4. The i ce-thickno s s parameter y\i expresses the nondirien- 
sional relation "between the ice thickness tj and thermal 
conductivity and the local heat- transfer coefficient h. 

If a certain size ice particle is assumed to have a nogligi- 
"blo effect on the propeller porformanco, the value of \i/ can 
"be determined at any point "by usin,^ the ice thickness and 
local value of h in equation (33), the value of the thermal 
conductivity of ice k.^ being a constant. The unit heat re- 
quired qi3g at any point at which \l/ is known can "be deter- 
mined hy the equation 

^bs = - % - ^A^ ^''^^ 



where is, in this case only, the blade-shoe surface 

temperature at the center of the smallest surface dimension 
(21) of the ice particle when the edge is at 32^ F. The 
use of a unidimcnsi onal analysis for a two-dimensional prob- 
lem is not considered to involve an appreciable error with 
the size ice particles and values of h normally encountered. 
The heat required at the blade-shoe surface is obtained by 
integrating the values of q-^g obtained by equation (34) as 

shown by equation (20), and the heat loss is determined in 
the same manner as was developed in the ice-prevention sec- 
tion of this report. 

A study of equations (32) and (33) indicates that when 
\l/ = 1 the ice-particle thickness and surface-temperature 
gradient are both zero. This indicates that when a v/ator 
drop strilves the blade shoe, it is immediately throv/n off as 
an ice flake. It is reasoned that there will be some mini- 
mum heat requirom.ent for propeller ice removal which will 
depend on the effect of an ice layer or ice particles on the 
propeller efficiency and it is possible that the optimum 
aerodyna,mic and thermal blade-shoe design will be for the 
condition of = 1. 
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APPLICATION 



It is apparent that in the application of the foregoing 
mathod to a specific propeller "blade-shoe design, the estima- 
tion of the effective angle of attack and resultant velocity 
at each pr opeller- Dl ade element is suhject to an error d\ie 
to body interference which will he reflocted in the calcula- 
tion of the heat-transfer coefficients. In addition, it is 
que st i onahl e whether it is judicious to design a propeller 
blade shoe for a single pr opell er - operat ing condition which 
fixes the blade angle and V/nD. With the foregoing in mind 
it is suggested that, in a specific design, the method pre- 
sented herein be repeated over a sufficient range of blade 
angles and values of V / nD to insure evaluation of the max- 
imum heat-- di st r ibut i on requirements. In this manner a fam^ 
ily of chordwise heat -d i st r ibut i on curves for each propeller- 
blade element is obtained over a safe propeller-operating 
range, from v/hich family envelope curves can be drawn for 
all elements to establish the required heat distribution over 
the entire blade shoe. 

In order to illustrate the method developed in this re- 
port, sample computations are presented hereinafter on the 
design of a blade shoe for a Hamilton Standard propeller 
blade l^o. 6477A-0. The computations presented are for one 
propeller-operating condition only and no allowance has been 
made for body interference. In an actual design, repetition 
of the method over a selected range of blade angles and val- 
ues of 'V'/nD, as indicated. in the preceding discussion, 
would be desirable. The Hamilton Standard blade No. 6477A-.0 
has been selected for the following analysis since this blade 
was used in the S-17F airplane propellers which were eauipped 
with blade shoes, tested, and reported by the NACA in refer- 
ence 2, The following flight conditions, for long-range 
cruising, will be used in the analysis: 



Ambient-air t em.peratur e , 0 

Pressure altitude, ft 10,000 

Indicated airspeed, mph , . . . . 155 

Propeller speed, rpm , 1012 



(Per these computations P is assumed to be at the 
V^/nD for maximum propeller efficiency.) 
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These are taken to "be the cr i t i cal c ond i t i ons for pro- 
peller ice protection "because the pr opell er-lDlade elements 
are near the maximum practical an,^le of attack. It is rea- 
soned that the high a ngl e- of-at t ack condition will result 
in the maximum heat requirement at the blade-shoe surface 
"because of the decreased extent of the laminar boundary lay- 
er on the forward face of the blade. 

The area covered by the blade shoe is assumed to extend 
chordwise to the 20-per cent-chord point on both surfaces and 
from the blade shank at the hub r/R = 0.145 to the station 
at v/hich the tip radius begins r/R = 0.942, v/here R = 69 
inches. The 20-percent chordwise coverage v/as satisfactory 
in the tests reported in reference 2, as was the radial 
bladc-shoc length which extended to r/R = 0.942. Actually 
the radial extent of the blade shoe would be determined by 
the amount of aerodynamic heating experienced by outer ele- 
ments of the blade. 

The computations follow the same sequence of steps as 
presented in the section Method. 

Determination of the Surface Heat -T ran sfe r Coefficients 

The chordv/ise heat -t ran sfer coefficients were determined 
at five blade stations (measured in inches from the center of 
rotation): 21.5, 25.5, 31.5, 42, and 63.5, which correspond 
to thickness ratios of 0.18, 0.15, 0.12, 0.09, and 0.06. 
The four outer-blade stations for the 6477A-0 blade were 
RAF-.6 airfoil sections, while the inner section was similar 
to a modified Clark Y airfoil section. 

Pr OT)eller-blade section velocity distribution .- The 
geometric angle of attack a was obtained by calculation of 
the helix angle $ and 'V'^/nD from the flight conditions. 
The values of ^-^^i "^^M ^^^^Q obtained from figure 10, 

reference 12. The blade angle 6 was obtained from the 
blade data sheet 157 of reference 12. The geometric blade 
angle 9 at 0.75R is given by the approximate relation 

Po.75 = ^0.75 + Gm - *M (35) 

By knowing the value of Po.75» relation of 0 to 

r/R was obtained by 
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P = 6 



^^0.75 ~ ^0,75^ 



(36) 



and a oy equation (l). The values of a, P> , and G 
are plotted in figure 5. 

The corrections for aspect ratio wore made by equations 
(2) and (3). The angle of zero lift and slope of the sec- 
tion lift curves were obtained from reference 13, and the 
section lift coefficient for each station was computed using 
e quat ion ( 4 ) . 

The velocity di st r ihut i on for the 18-per cent-thi ck sec- 
tion was ohtained hy the use of the method presented in ref- 
erence 5, An IIACA 0015 airfoil was used as the reference 
profile hecause the leading-edge radius was similar to that 
of the modified section. The velocity distribution over the 
NACA 0015 airfoil is well defined "by experiment. The values 
of "^r /""^R v;ero obtained as indicated in reference 5 by 



(37) 



The velocity distri- 
velocity "^^r/^^R* 



where P^, is the of reference 14. 

butions vrere obta.ined by adding to this 

change in velocity AV/Vp^ due to the difference in shape 

between the reference profile and the 18~per cent -thi ck sec- 
tion 



then 



AV 

V 



P. 



(38) 







F/4 




^ ^R 


^ ^f/^R 




_ 


P/4 


Vr 




^f/v^R 



(39) 



where P is the pressure corresponding to the normal force 
at a chord point given by the equation 
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(40) 



as defined in reference 5. The velocity di st r ilDut i on over 
the 12-per cent-thi ck RAF-G blade section was obtained by the 
method of reference 4 and was modified for the various val- 
ues of thickness and camber to give the velocity distribu- 
tions for the remainder of the blade sections. Veloc- 
ity dif trubutions for the five blade stations are plotted in 
figures 6 to 10. 

Boundary>-layer heat transfer .- The chordwise position 
of the laminar- separ at i on point v/as determined by approxima^*. 
ing the actual velocity distribution by a double-roof profile 
and finding ^g/'^M ^'^oia figure 2. The approximate double- 
roof velocity profiles are shown as dashed lines for the 
upper and lovjer surfaces on the velocity-distribution curves 
in figures 6 to 10. The double-roof profiles are drawn so 
as to give approximately the same lift (area beneath the 
curve) up to the laminar- separ at i on point as the velocity- 
distribution curves. Considerable judgment must be used in 
the dra^^^ing of the profiles in order to obtain reasonable 
results. Theory indicates that, for the velocity distribu- 
tion shov/n in figure 10, laminar separation v;ill not occur 
on the lower surface of the propeller-blade element. Exper- 
ience, however, indicates that laminar separation is likely 
to occur and, since the assumption that separation does occur 
is conservative, the laminar- separ at i on point shown for the 
lower surface in figure 10 has been used in this analysis. 

The values of boundary-layer thickness and heat-trans- 
fer coefficients v/ere computed by equations (6) to (14) and 
are plotted in figures 6 to 10 with the corresponding veloc- 
ity distributions. Because of the changes in lift coeffi- 
cient of the blade sections as the blade angle is changed, 
the location at which the heat transfer (due to forced con- 
vection) approaches zero will move fore and aft from the 
design position. In order to protect the region of the 
blade at v/hich lev; values of h occur at other than design 
Cj, the curves of h as related to x/c were faired as 

shown in figure 11 and the faired values v/era used in suc- 
ceeding coraputat ions . 
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Heat Required for Tee Prevention 

H eat req uired to maintai n the "blade- shoe surfa ce at 

3 2 F.- Equations (18) and (19) were used to determine the 

unit pov/er input to the "blade shoes. The water content of 
the air stream v/as taken to "be 1.84 grams per cu"bic meter, 
tho maximum water content moasured in clouds hy Kohler, as 
descri"bed in reference 15. 

When the various factors in equations (18) and (18) are 
fixed by the design conditions, the equations "becon-e 



^^ts U,. = ^ta^ - 0.833 X 10-^ V^/-] 



'ystag 



+ 13.1 - 8.18 X 10"^' Vjj^ 



(41) 



and 



..-4 



q,^^ = h[32 - 0.667 X 10 ] 



+ Hi [13.1 . 8.18 X 10 ^ Vp/] (42) 



The heo.t per unit area required to X'i''0^''ide ice protec- 
tion on the surface of the "blade shoe q^g can then "be eval- 
uated "by suhsti tut ing in eqiiations (41) and (42) the proper 



values of h (from fige ll), , and Vp_ for the respective 
point on the "blade-shoe surface under consideration. The 
total heat req^uired at the hlade-shoe surface to provide ice 
protection is obtained "by integrating equation (20) over the 
hlade-shoe area, from r/R - 0.145 to r/R - 0.942 and from 
^4. n X. x/c - 0.,2. For the three propeller "blade 



at x/c-0 to 



^bs 



0.942 



^. 3R 



0 0 14 5 



X / C =0.2 



r 



ds 



(43) 
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The heat req_nired at the 'bl-'ide-shoe surface to provide 
ice protection v;as ohtained "by pcrforning the integration oi 
equation (45) "by ^;raphical Lieans, The power required to pro 
vent ice formations on a Hamilton Standard propeller v/ith 
three lTo» 6477A-0 blades expressed in electrical units is 
3957 watts at the desi.f:n conditions v/ith the losses through 
the after portion of the "blades neglected. The results of 
the computations of the unit power distribution are ^^'iven in 
f i gu r e 1 2 • 

Heat loss r- The heat loss from the aft portion of the 
"blade is assumed to occur over the rear 80-^percent chord, 
and from the blade root r/H =: 0^145 to the tip region r/R 
0.942. With these limits, eqiiation (28) becomes 



0.942 



0.14 5 



r' 



1 c 0 



.7 



■^as 



)tt dx 



0 c 2 



1 o 0 



0.2 



(44) 



The thickness of the orter layer of the hlade shoe for this 
analysis is 0.020 inch or 0. 00167 foot, v/hile the inner 
layer tj is 0,040 inch or 0,0033 foot thick. At any hlade 
station, the values of S^^/S^ and S^^/S^ are taken to he 

0.4 and 1,5, respectively (fi^. 1 ) ; also, from the wedfje 
hlade element 



^''^jrax L 



(45) 



where h is measured fron the hase of the wedf^e and (t/c) 
is the maximum thickness of the "blade section. Therefore,' 
equation (27) "becomes 

0.0133 + 4Ati 



- T, 



Ta) = qas 



t/c 



1.6h 



1.25 b/c] 



n 



(45) 
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and v;hen 



and 



kj = 0.08 Btu/hr, ft^, °F/ft 
k-^ = 117 Btu/hr, ft^, °F/ft 



then, inserting the foregoinfV values of thermal conductivity 
and rearranging the terms, equation (45) 'becomes 



To - 



'•as 



0.167 + 50At, + 7-7-T 2.^137h + 1 

(t/c)^^ [1 - 1«25 (d/c)] h 



(47) 



3y siiccessive approximations, values of Atj, Q,g ^ , and 
'■^"bs which are in reasonable agreement can "be found and the 
heat loss through the after portion of the propeller can he 
evaluated "by integrating equa,tion (44)^ A graphical inte- 
gration of this equation indicates a total heat loss for the 
entire propeller of 1109 watts* 

The total po\/er required to prevent ice formations on a 
propeller v/ith three 6477A-0 hlades v/ith the di 3 1 rilDut i on 
shov/n in figure 12 is 395 7 + 1109 = 5 06 6 v;attSo 

The calctilated povrer, 5066 v/atts, .is greater than the 
pov/er used diiring the flight tests in natural icing condi- 
tions reported in reference 2 and is greater than the avail- 
able pov/er from generators which are now under development 
for propeller ice prevention^ Although the process of ice 
prevention is probably the ultimate objective in the develop- 
ment of protection for propellers, the process of ice removal 
as indicated by the results of roforencG 2 appears to be a 
practical solution more suitable for immediate applicatione 



Heat Required for Ice PweriOval 

The heat required for ice removal will depend on the 
maximum size of ice particle that is allowed to form on the 
propeller blade. Equation (32), when T^ = 0^ F and the 

factor [l/--cl^cl^ evaluated for the particular blade- 

shoe design, becom.es 
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1 .iL^jir^^^' (48) 

and since the heat required to remove an ice particle v;ill 
"be a function of the temperature gradient at its edge, the 
values of T^/T^^ for the suhject propeller and blade shoe 
were coiiiputed and plotted in fi.-ure 13 for various values of 
\|/ and h"^l. 



In order to clarify the effect of the surface tempera-- 
ture, an ice particle 0^3 inch in diameter and of variable 
thickness v/as assumed to "be attached to various points on 
the "blade shoe* The chan.^e in air flow over the "blade duo 
to the ico particle is assumed to have a negligihle effect 
on the heat transfer^ The heat required was "based on T. 

2 * °° 

for the various values of \|/ . 

The curves of figure 14 were computed to shov; the heat 
required to remove the assumed ice particle as a function of 
thickness and heat-transfer coefficient. In order to oh tain 
the correct value for the unit pov;er at an amhicnt-air tem- 
perature of 0^ at any point on tlie hlade, the unit power 
from figure 14 must he corrected oy a function of the aero- 
dynamic heating as shown "by the enuation 

With an3^ ice particle of smaller diameter than that as- 
sumed, more heat or a greater thickness of ice v;ould he re- 
quired for its removals Observations during icing flights 
have shown that in icing conditions at low temperatures 
(O to 10^ ?), ice forms in a very narrow ridge, sometimes 
as narrov; as l/4 inch^ along the hlade leading edge. This 
condition has heen substantiated oy other observers (refer- 
ence l). The ice which forms at these low temperatures is 
particularly hard and tenacious* In viev/ of tlEic occurrence 
of such icing conditions, the use of less heat than that re- 
quired to maintain the blade-shoe surface at 32'^ ? ( = l) 
in dry air, seems to be inadvisable. 

The distribution of the unit pov/er to obtain ico pro- 
tection by removal at all points on the blade-shoe surface 
simultaneously, based on figure 14 at = 1, is plotted in 
figure 15, Figure 16 shows the required unit heat distri- 
bu.tion on the developed blade-shoo surface. 
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Graphical integration of the unit heat required over the 
"blade-shoe area gives a value of 1750 watts required to effect 
ice removal at \|/ = 1 . The heat loss through the after por- 
tion of the iDlade was computed through the use of equations 
(22), (23), (43), and (44) and was found to "be 640 watts. The 
total power required for ice protection by the removal process 
at ^1/ = 1 is therefore 2390 watts. 



DISCUSSION 

Benefits Obtained from the 0ptim>im Heat Distribution 



In view of the relative simplicity of construction of a 
propeller-blade shoe with a stepped type of heat distribution, 
such as was used in the tests reported in reference 2, when 
compared with a blade shoe having an optimum heat distribu- 
tion as shovm in figure 16, it is important that the blade- 
shoe designer realize the benefits which may be expected by 
utilizing a blade shoe having an optimum, heat distribution. 
To indicate the advantages of an optimum heat distribution, 
the method of analysis developed in this report has been 
employed to compare the amount of power required by the two 
types of distribution to provide the same degree of ice pro- 
tection and to determ.ino the thickness of ice which would be 
built -up upon a blade shoe with a stepped heat distribution 
utilizing a power input equal to that required for complete 
ice removal on a blade shoe with an optimum heat distribu- 
tion. In this comparison the optimum heat distribution 
shown in figure 16 has been compared with a stepped heat 
distribution in which the distribution over the leading-edge 
third of the blade shoe is double that over the remaining 
area. The chordwise heat distribution for the blade shoe 
with tho stopped distribution has boon considered constant 
over the radial extent of the blade shoe and has been based 
upon the point of maximum heat requirement for the blade- 
shoe surface. The chordv/ise distribution for each type of 
heat distribution is shown in figure 17 at the critical pro- 
peller station. 

The calculated total heat required for the distribution 
shown in figure 16 at T^ = 0^ F and \l/ = 1 is 2390 watts, 

while the calculated power required for the stepped distribu- 
tion, based on tho point of maximum heat requirement on tho 
blade-shoe, surface (station 42 stagnation region) at the 
same conditions, is approximately 4800 watts. To enable a 
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comparison i/ith the e:cperiment al results, the heat renuired 
"by the optihaim di s tr i"biiti on and the stepped distrihuti on v;or 
computed for the test conditions^ These conditions v/ero con 
paralDlo to those listed under the section of this report des 
ignated Application, except for the ambient-air temperature 
v/hich v/as 10^ F instead of 0^ The computations indicate 

that 1C40 v;atts are required by the optimum distribution and 
approximately 3500 watts by the stepped distribution (both 
at \|/ = 1 ) , The use of an optimum heat distribution rather 
than the stepped distribution should, therefore, reduce the 
total heat required by approximately one--half for the same 
decree of protection af^?:ainst ice formation. Lacking more ex 
tensive experimental confirmation, tho absolute magnitude of 
the foregoinj.'; po\\rer requirements are r>ubject to some doubt; 
hov/ever, their relative ma^'-^ini tude should be reliable* 

If the stepped heat distribution is used, but with a 
power input of 2390 watts at tho flight conditions listed, 
ice v/ill accumulate until its thickness supplies sufficient 
insulation for the removal process to take place. Applica- 
tion of the foregoing method to compute this thickness indi- 
cates that an ice cap 3c 42 inches hi.;^h by 0o3 inch v;ide 
(width assumed, based on observations) would form before 
being removed by melting at the base. Obviously, such an 
ice structure would be unstable and break, leaving the lead- 
ing edge v/ith a rough, broken ice cap and, thus, the blade 
shoe v;ould fail in its function of ice removal,. 



Comparison with Experiment 

The foregoing deductions r.ust be qualified and consid- 
eration must be given to the lack of experimental verifica- 
tion of the proposed methods Such experimental verifica- 
tion is particularly needed because of the assumptions v/hich 
were necessary in the development of the computation proce- 
dure. The results of reference 2 are the only experimental 
data that can be compared with the analytical results. 
These tests indicated that satisfactory ice-removal charac- 
teristics would be obtained with 2100 watts applied to tho 
stopped distribution of the blade shoes tested not known), 
while the computations indicate that approximately 3300 v/atts 
v/ill be required with the same distribution at similar con- 
ditions (v = l). In view of this lack of agreement, a crit- 
ical view of the assumptions upon v/hich the method is based 
is in order o 
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Review of Assumptions 

As the riethod v;as applied, two ac^ sumpt i ens appear to oe 
conservative. In the first place, "body interference was ne{,'^:- 
lectedj which eliminates from consideration the "blocking of 
the air and the resultant reduced velocity over the region of 
the propeller close to the shanks Data are available which 
indicate that, on an air-cooled nacelle-propeller comlDination 
of the same general dimensions as that on the B-l?? airplane, 
the axial velocity is approximately free stream from the tip 
to the 0o5 r/R station and then decreases linearly to ap- 
proximately one--half free--sbream velocity at the 0^2 r/R 
station* A reduction in the comT)uted heat requirements for 
the inner-Dlade sections would ohviousiy occur if allowance 
v/ere made for the reduced velocity* 

A secondary result of the reduced axial velocity will 
he an increo.sed effective olade angle over the inner stations, 
'The magnitude or direction of the deviation req\xlting from 
this difference is difficult to predict, "but these "blade sec- 
tions will Cc?.rry a higher lift coefficient than that assumed. 
By reference to figures 5, 7, and 8, it can he reasoned that 
this will Cause an increased peak velocity on the upper sur- 
face, thus decreasing the extent of laminar flov/, hut v/ill 
have a compensating effect on the lower surface v/here the 
peak will tend to he lowered* 'A^aich of these changes will 
predominate could only he determined hy a detailed calcula- 
tion. 

A second conservative assumption is the neglect of cen- 
trifugal-force effects. This was justified originally on 
the hasis of the possioility of a reasonably uniform glaze- 
ice layer formed from the maximum possible water content of 
the air stream at an ambient-air temperature of 0^ F. Under 
these conditions I the ad-hesion of ice to the blade-shoe mate- 
rial is greater than the strength of the ice itself. That 
this assumption may be unduly rigorous can be deduced from 
the result of reference 2, which indicated the possibility 
that no heat at all might be necessary to de-ice the outer 
portion of the blade. VHiether this is due to aerodynamic 
heating or centrifugal-force effects is difficult to deter- 
mine. Kov/ever, the basis for the computations of the aero- 
dynamic heating is reasonably sound, while the possibility 
definitely exists that an ice formation more porous than the 
glaze ice assumed v;ould adhere less strongly and, therefore, 
be thrown off by centrifugal force; also, the porosity of 
this ice layer should have a marked effect on its thermal 
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resistance. ITo answer to these possi"bili ties can "be proposed 
without further study. 

Several unconservative assumptions were made in the anal- 
ysis which would tend to cancel the error caused "by the ne- 
f^lect of "body interference and centrifugal force. The r.ore 
important of these is the assumption that transition v/ill 
occur "by the process of separation in the adverse pressure 
gradients Actually if any ice forms around the leadin,^ edge, 
it is likely to cairse local peak pressures or, through roug!?.- 
ness, promote earlier transitionso 'Thus, instead of the 
fairly extensive lamina.r flov; over the upper surface, tiirhu- 
lent flov; might "be present. Since the acciiracy of the method 
is largely dependent on the computed location of "boundary- 
layer transition, this would have an important influence 
tending to make the method unconservative o Tor ice accumula- 
tions occurring with ice removal at if of less than 1, the 
actual formation of ice layers and ice flecks will vary 
widely over the "blade and will tend to increase turbulence 
further and to effect the rate of heat transfer. 

The effects of humidity^ evaporation, and degree of 
supercooling of the v/ater drops have also "been neglected; 
hov/ever, these effects are considered to he sufficiently 
small as to he compensated "by other conservative assumptions. 



Reconmendat i ons 

The c[uantity of pov;er required for ice prevention (5055 
watts) appears to "bo prohibitive at the present stage of de- 
velopment. Because of the weight of the huh generator that 
would oe required, the advantage of the lighter v/eight of 
the th prmpj. --elect ri c system in comparison with other methods 
would he lost* The m.05t expedient present solution appears 
to "bo a design for ice rem.oval at \\f equal to lo If advan- 
tage is taken of the reduction in heat requirements made pos- 
sible by use of an optimum heat distribution, and if the 
thickness of the conducting layer is kept to a minimum to 
rediice the heat lost through the after surface of the blade, 
a blade shoe which will provide satisfactory ice protection 
with available hub generators should result.. In order to ob- 
tain the required heat distribution with a satisfactorily 
thin conducting layer, a m-aterial of varying electrical re- 
sistancoj but constant thickness, v/ould be advantageous. 

As a resiilt of limited tests, it has been found that 
present conducting materials do not have satisfactory 
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resistance to a"brasion and wear. It haf? "been thought that a 
thin layer of nonconducting ru*bher over the entire hlade- 
shoe surface would provide the al^rasion resistance necessary 
for satisfactory service life. The effect of the addition 
of a protective layer can "be determined hy including its 
thickness in the thickness of the conducting layerr, 
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APPENDIX 



An expression for the temperature gradient in the edge 
of an ice layer is to he developed for determining the size 
and thickness of an ice jparticle that v/ill adhere to a blade- 
shoe surface v/ith different \.init power inputs. In the deri- 
vation of this expression symhols will "be used which do not 
appear in the nomenclatiire of this report (C^^, C^, f, vi, 
and v)o It v;as felt that hotter continuity and clarity 
could be had by defining the symbols as ^hey appear g 

The diagrams of figure 4 illustrate a section through 
the propeller blade normal to the edge of an ice layer and 
the temporaturo proiilco 
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From 1=0 to I = the extent of the ice layer, 

then Gi:u:plifying and dividing "by dl gives 



dCij dq^^^ 



dQ. 



dl 



d I 



dl 



Furthermore 



n " "crc I . 7 



and differentiating £;ivos 



2 

d T- 



dl " ^"cl^cl :2 



When the pov/er input is considered to supply Q,^^ then 

dC = 0" dl 



and from figure 4 



m- 



1+ ^ 



! dl 



Substitution »sives the equation 



2 

d T, 



h ic^ 



Dividing through "by (k^jt^j) the equation takes the forn 



d^Ti 



(Ti - To) - ^ 



where 
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If 



then 



u = 



hk, 



2 

a I 



= U f - V 



The solution of this eqiicttion is of the form 



f - C-e + C^e + 



where 0;j_ and are cons t ant of integration and 



^i ^^^i 
hkl 



00 0 



Suh stitution gives 



and when I "becomes large (T-j^ - ?q) approaches a finite 
value asymptotically as indicated in figure 4» Therefore, 



and 



= 0 



Ti - Cse + 



When I approaches zero 

T. 
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Similarly, the equation for the surface temperature from 
1=0 to -00 takes the form 



and 



d 

d I 



^1 = C.e^^ + 



-00 



Wiier e 



u 



h 1 



At 1 = 0 



and 



Also, v;h.en the sulDscriptr. + and - are introd-'aced to indi- 
cate positive and negative values of I, 



d? 



dl 



and 



11 I 



When I eq^iials zero. 



dT. 



dl 



and 
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Squatin^^ the two foregoing equations and siil) s t i tut ing the 
value for C^^ result in the expression 

-u^Cs = u^(Cs + - ) 

and, r earran^'^ing terms. 



' CO 



U 



^- 1 



Substitution of the equation for in the equation for 

T^^ , for 1 = 0 to I - CD^ results in the expression 



1. 



— + 1 



Because o-^^ is uniform ovor tho element of "blade surface 



being considei-ed, 



or 



) = (T, - T,) 

k. + ht. " 



k. + lit, ^ 



: ) + T 

0 ' o 



and, making the final sulistitutions in the exprcjssion for 



k. 



1 (T - T. ) + T 

k, + ht. 0= 0 



^i-^i-^^i '-'-cI^cJJ 



k4 + ht. 



L-^i 



+ 1 



Introducing the ice-thickn3ss parameter 



^k. + ht. ^ 
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the foregoing equation may "be written in the form 




and when T = 0 

0 



T. 



= 1 - 



L - 



r I 



1 2 



P 1 




Figure 2.- The ratio of the local velocity at laminar separation, V-t to the mflxi""^"' local velocity t 
Vjj, as a function of the ratio of the slopes of the double- roof velocity profile, p« 




Figure 3.- Variation of f[^]with ^ for the turbulent boundary- 
layer thiclaiefls calculations. 




Figure 4.- A section view and temperature profile normal 
to the edge of an ice layer. 
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Figure 5.- Variation of propeller— blade angle, B. helix angle, $ , blade- 
angle distribution, 6, and angle of attack, a, with percent 
radius for the Hamilton Standard blade No. 6477A-0 at the design conditions. 
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FIGURE 12, - UNIT POWER QlSTRlBUTlOM FOR COMPLETE 
ICE PREVEMTION. 
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Heat required to remove an ice particle .3 inch in diameter ""^ 
(Tq = 0° F) , neglecting aerodynamic heating. 
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Figure 17.- Comparison of the stepped and optimom unit power 
distributions at station 42. 



